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Abstract 

In this paper we emphasize some conceptual points related to the kinetic foundations of 
relativistic hydrodynamics. We summarize previous work and focus on the construction of 
the heat flux from a kinetic theory point of view. A thorough discussion addressing aspects 
concerning stability, causality and the construction of an appropriate stress-energy tensor 
is included. 
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I. INTRODUCTION 



The first solid theory formulated to deal with irreversible processes is due to Lars 
Onsager jl| and its first version was published almost eighty years ago. Extensions 
and modifications of his method were made by Casimir |3j and almost simultaneously 

n 

by Meixner |4j, all of which gave rise to what we presently know as linear irreversible 
thermodynamics (LIT). However, in 1940, before LIT was finally accomplished, C. 
Eckart published three papers on irreversible thermodynamics [5] , the first two more 
or less along the lines of LIT and the third one dealing with the relativistic irreversible 
thermodynamics of a simple fluid. Thus arises the question: is this work the correct 
extension of LIT to relativistic systems? The answer has been subject to continuous 
debate for over forty years and this is precisely the motivation of this work. It is 
important to point put that no new results are included in this publication. Following 
a very much appreciated suggestion from an anonymus referee, this review article has 
been written as a thorough discussion of some key theoretical aspects that lie beneath 
recent results that can be found in the literature. Some basic equations are included 
for the sake of clarity. 

It has recently been shown [6] that the source of the generic instabilities of a 
relativistic fluid |7| is the misuse of the constitutive equation for the heat flow first 



proposed by Eckart in 1940 [5]. In that work it is clearly shown how the acceleration 
term in the expression 

Jh = ~*K { T '" + (!) 

is the one that leads to the exponential growth of fluctuations around the equilibrium 
state in the linearized system of relativistic transport equations. Here JYq, is the heat 
flux, k a relativistic thermal conductivity, h^ v is the usual spatial projector defined 
as hF" = g^ u + for a (+ + H — ) signature, T is the local temperature, c the speed 



of light and the hydrodynamic four velocity. The covariant derivative is denoted 
by a semicolon, the ordinary one as a comma and a dot implies a total proper-time 
derivative such that a u = u^a v ^. In this work we point out two major drawbacks of 
this constitutive equation and thus proceed to address the matter from the kinetic 
theory point of view. We repeat some key equations leading to the constitutive 
relation obtained from the relativistic Boltzmann equation as well as the ones needed 
in order to address both the violation of the linear regression assumption and the 
causality of the first order (in the gradients) theory. 

To accomplish this task in Section II we start by recalling the basic elements of 
Onsager's theory in order to clearly expose the shortcomings introduced in first order 
theories by Eq. ([TJ). In Section III we outline the calculation of Ref. |8j where a 
constitutive equation for the heat flux is derived by solving the Boltzmann kinetic 
equation for a relativistic ideal gas in a BGK-like approximation. In Section IV we 
briefly summarize the key steps and arguments favoring the kinetic theory result 
and showing its consistency with linear irreversible thermodynamics. Finally, the 
key aspects of the work are summarized and discussed in Section V. 

II. GENERIC INSTABILITIES AND ONSAGER'S HYPOTHESIS 

The constitutive equation (CD), arises from phenomenologically enforcing the sec- 
ond law of thermodynamics to a system in which the stress-energy tensor includes 
relativistic corrections terms proportional to the heat flux. Such tensor is obtained 
as the most general decomposition of a second rank tensor in this framework. This 
procedure is the standard one in non-equilibrium thermodynamics, however it leads 
to a consitutive equation which is at odds with its basic tenets in two senses. First, 
according to this theory, fluxes must be coupled to the forces given by the gradi- 
ents of the intensive variables. This requirement is clearly not met by the second 
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term in Eq.flT]), the acceleration is neither an independent variable nor a thermo- 
dynamic force. The second problem, as mentioned above, is the fact that it leads 
to a violation of the linear regression assumption, a behavior usually referred to as 
"generic intability". As was pointed out in Ref. (6|, Onsager's hypothesis states that 
spontaneous fluctuations of the state variables around the equilibrium state, whose 
origin is purely microscopic, should relax following the linearized equations for such 
variables. Thus, the exponential growth obtained by Hiscock and Limdblom is really 
a violation of this assumption which is key in the construction of irreversible thermo- 
dynamics and in particular in the proof of Onsager's reciprocity relations. Indeed, 
the theory of irreversible processes developed by Onsager seeks the establishment of 
a connection between these processes and the spontaneous fluctuations that due to 
their microscopic nature appear in the thermodynamic variables of equilibrium sys- 
tems. To do so he used the local equilibrium assumption (LEA) \^ flo| by assuming 
the thermodynamic state of a system to be determined through a set of extensive 
variables (a\, a.21 a n ) = a (energy, volume, particle number, etc) with which all 
thermodynamic variables, in particular the entropy S, may be defined. These a's are 
further redefined so that they vanish in the equilibrium state, 5* {p) Eq = S ((Pj = Sq 
and thus, the thermodynamic forces acting on the system to seek equilibrium are: 



while the fluxes are simply the time derivatives of the a's. The essential assumption 
in Onsager's formulation is that the fluxes are linearly related to the forces Aj so 
that 



Further, to determine the force Xj it is assumed that S (a) can be expanded in 
a Taylor series around equilibrium. Since the linear term vanishes and defining 




(2) 



ddi 
~dt 




n 



(3) 
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AS = S (a) — So one gets that 

j n n 

AS = -^EE^«»«j ( 4 ) 
z i=i i=i 

where the g^s are equilibrium properties. Using Eqs. (121) and (J4]) it follows immedi- 
ately that the linear laws take the form: 

Xj 1^*3 + 9ij a jj =0 i = l,..,n (5) 

where Rij = (Lij)^ 1 . Assuming that the spontaneous fluctuations around the equilib- 
rium state of any system obey the same equations as the macroscopic variables, e.g. 
Eq. (j3J), using the methods of statistical mechanics and the fact that the microscopic 
equations of motion are invariant under time inversions (microscopic reversibility), 
Onsager was able to show that the matrix L^, the transport coefficients matrix, is 
symmetric 

Lij — Lji (6) 



The proof o 

n 

chanics [11|- 



' this statement may be found in any good textbook in statistical me- 
14} . or in Refs. What is of outmost importance is to understand 

that once we admit the validity of Onsager's Reciprocity Theorem, Eq. (jSJ), we 



are implicitly admitting the linear regression of the fluctuations hypothesis, which 
has been thoroughly corroborated by experiments in colloid systems [16], light 
scattering of fluids [171 and other systems Jl8|. 

However, Casimir 3fl and J. Meixner [4] independently pointed out that many 
irreversible processes did not fit into Onsager's formulation, in particular the cases 
where the state variables are field variables, such as in heat and electrical conduction. 
In such cases the state variables satisfy conservation equations and the forces are 
given by gradients of intensive variables instead of defined as in Eq. (j2j). Two major 
changes appear: Eq. (JI} is substituted by the entropy balance equation 

p— + divJ s = a (7) 
at 



which results strictly from the validity of the LEA and the conservation equations. 
In Eq. ([7j) p is the mass density, o is the entropy production and J s is the entropy 
flux, which is a measure of the entropy flowing through the boundaries of the system. 
Secondly, for an isotropic system a turns out to be a sum of products of forces and 
fluxes of the same tensorial rank. Thus, the linear relationship in Eq. is consistent 
with a > but now the fluxes are no longer the time derivatives of the field variables 
but the currents. Physically, a is just the local version of what Clausius in his last 
paper on thermodynamics called the "uncompensated heat", a measure of the heat 

BLerated by the dissipative effects occurring in a process taking place in the system 
|. Following the linear relation hypothesis, the entropy production turns out to 
be a quadratic form in terms of the forces, usually written as 

n n 

a = J2T, L ikXiX k (8) 

i=l k=l 

such that, if the transport coefficients Lik are positive, then a > 0. 

Onsager's reciprocity theorem, stated in Eq. (Q, is kept as a potulate in the 
extensions of his theory that constitute nowadays LIT formalisms. This means that 
one is ultimately also accepting the linear regression of fluctuations hypothesis. To 
summarize this section, our objection to the use of Eq. ([I]) comes from the fact that 
it is inconsistent with classical irreversible thermodynamics in two senses: it does not 
have the structure given by Eq. ([8]) and it leads to a violation of the linear regression 
of fluctuations hypothesis. The proof of this last statement follows directly from the 
fluctuaction analysis using standard hydrodynamics techniques. The full calculation 
is shown in Ref. |6fl. 
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III. THE KINETIC THEORY APPROACH 



The contradictions that arise from coupling the hydro dymamic acceleration, or 
any type of term that is not a gradient of a state variable, with the heat flux have 
been clearly stated in the previous section. Since this coupling seems to arise from 
following a phenomeno logical treatment, one is naturally lead to the question of 
whether such a relation is predicted by kinetic theory which, as in the classical case 
provides a microscopic framework for the establishment of constitutive equations 
from basic principles. The procedure is indeed well-known in the non-relativistic case 
through the use of Boltzmann's equation for a dilute gas. In the relativistic 
long as the detailed expressions for the transport coefficients are not required, one 
can accomplish this task by using the simplified form of the relativistic Boltzmann 
equation known as Marie's equation 20]. This equation reads, 



f _ f (o) 

v a f, a = - J -^— (9) 

T 

where v a is the molecular four- velocity, / is the single particle distribution func- 
tion and /(°) the equilibrium distribution, i. e. the solution to the homogeneous 
Boltzmann equation. In Eq. (J9]) r _1 is a parameter which contains the information 
arising form the collisions between the particles. However, as mentioned above, this 
information is not required for our purposes. Such details are required in order to 
calculate transport coefficientes while the structure, in the sense of the dependence 
with the state variables, of the dissipative fluxes can be assesed without them. 
The solution to Eq. ([9]) can be obtained by the standard methods of kinetic theory, 



which have been shown to be valid also in the relativistic case in Ref. |2l| . Here we 
only outline the procedure since the full calculation has been shown elsewhere 
In order to obtain the solution to first order in the gradients one uses the Chapman- 
Enskog method in which the distribution function is proposed as 



7 



/ = / (0) (1 + 



(10) 



The first term corresponds to the Euler relativistic regime and the second one, where 
<p is the first order correction in the Knudsen parameter, to the Navier-Stokes level. 



The equilibrium distribution function in Eqs. (jl 
which, in the non-degenerate case, is given by 



and 
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ITU]) is a relativistic Maxwellian 



22J 



/ 



(o) 



?? 



-e zc 2 



11) 



4vrc 3 2/C 2 (i) 

where m is the rest mass of the particles, z = the usual relativistic parameter 
and K, n the modified Bessel function of the n-th kind. Subtitution of this hypothesis 
in Eq. ([9]) leads to 



-TV 



d f(0) df (fl) 



-n.a + 



d f(o) 



duP~ U > a , 



(12) 



dn ' dT 

The derivatives can be readily calulated and subtituted. The Euler (previous or- 
der) equations are used to write the time derivatives, appearing in the sum when 
a = 4, in terms of the spatial gradients. Notice that this procedure will naturally 
lead to a coupling of the dissipative fluxes, calculated as moments of 0, with the 
thermodynamic forces. The solution, in the fluid's comoving frame can be written 
as 



/ 



(0) 



'TV f 



cp zc 1 




2zK 2 (I) 2zK 2 (i) 



(13) 
(14) 



where the isothermal compressibility, /3 the thermal expansion coefficient and 
p = (ne + p) /c 2 where e and p are the internal energy and pressure respectively 
which, in turn, have to still be expressed in terms of the scalar indepentend variables, 



8 



T and n, through some equation of state. This is required for consistency, since the 
chosen representation is n, T and u u . Then, the forces present in the solution <fi are 
exclusivelvthe temperature and density gradients. By means of the Enskog transport 
equation [8], the heat flux in the comoving frame can be shown to be 



n (mc 2 v^ = mc 2 / v^f^^jdv 



J [Q] 



such that, since f^(p is coupled with VT and Vn, the corresponding constitutive 
equation has the following form 



[Q] 



-K (L TT T'» + L nT nn 



(15) 



where L?t is an "effective relativistic thermal conductivity". The new transport 
coefficient L n T has no classical counterpart. We would like to remark that e quat ions 



similar in structure to Eq. (fT5|) were proposed by previous authors 20|, |22|, |23| but 
also to strongly emphasize that the term proportional to the gradient of n arises 
here since the pressure p is not an independent variable, and thus the rules of linear 
irreversible thermodynamics compels us to express Vp in terms of Vn and VT. 

The fact that kinetic theory predicts a heat flux that depends only on gradients of 
state variables is widely accepted. However, as mentioned above, the phenomenology 
seems to relate it with the acceleration. This ambiguety remains an open issue that 



will be addressed elsewhere. In the rest of t 
results given by Eq. (fT5l) or variations of it 



lis work we focus on the kinetic theory 
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IV. RELATIVISTIC IRREVERSIBLE THERMODYNAMICS: CAUSALITY 
AND STABILITY 

As a natural extension of the arguments summarized in Sect. II, it is intuitively 
clear that the formulation of relativistic irreversible thermodynamics should start 
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by considering systems whose states are described through field variables which are 
continuous functions of the space-time coordinates x a = (x±, X2,x^, ct) and will hence 
satisfy conservation equations. If N a = nu a is the particle flux, 

N% = (16) 

is the statement of conservation of particles, i. e. the continuity equation. For 
the other two state variables, u v and e (or T) the balance equations are no longer 
independent. Moreover, the resulting equations must imply the local relativistic 
versions of the first and second laws of thermodynamics. The energy-momentum 
conservation law 

= (17) 

can be obtained directly from Boltzmann's equation. However, the still open question 
is whether this is to be matched with the conservation law implyied by Einstein's 
field equation 

Rfjiv — -^Rg^u = kT^u (18) 

where R^ — \Rg llu is Einstein's tensor describing the geometry of space-time. In 
Eq. ( TT8l) . T^u is the energy-momentum tensor which accounts for the properties of 
matter and k is the coupling constant. Since the covariant derivative of the Einstein 
tensor vanishes one has a conservation law T.$ v = 0. As already shown by Einstein 



himself 



24j . the Euler equations of relativistic hydrodynamics for an inviscid fluid 



follow directly when 



TIF 

= —u» u v + ph^ (19) 



However, he never address the dissipative case. 

On the other hand, from a purely hydrodynamic point of view and in order to 
include heat in the total energy account, Eckart resorted to an irreducible decom- 
position for T^ u and identified relativistic corrective terms proportional to the heat 
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flux. His proposal, in his own words, was that no assumptions would be made for 
the form of T^ u , but it would be used to define other quantities such as the internal 
energy e and the heat flux Jjgj. He thus obtained 

TIF 1 

T*" = + ph^ + W v + ^(Jfou" + Jfou") (20) 

although he never mentioned whether Eq. (I2UI) should be the one compatible with 
Eq. (|18|) . One should be cautious in undertaking this step since heat is not a state 
variable, but energy in transit. An alternative to this proposal was examined in 
Ref. [25] following Meixner's ideas by omitting the last two terms in Eq. (1201) and 
introducing heat in the definition of a total energy flux in a similar fashion as done in 
non-relativistic LIT. This approach led to a debate regarding the question of whether 



the heat s 
tensor 



26] 



lould be included in Einstein's field equation through the stress-energy 



23, • 



In this work as in Ref. |8fl and Section IV of Ref. 28|, the tensor proposed by 
Eckart is used while Meixner's formalism is kept as a separate alternative. Indeed, 
introducing the constitutive equation ( 1T5I) in the set obtained from ( 1T6]) . (fl7|) and 
([211]) a closed set of equations is obtained which will not be repeated here. The next 
step in order to examine if Onsager's linear regression of fluctuations hypothesis 
holds true implies linearizing such system by setting T = Tq + 5T, n = no + 5n and 
u k = 5u k (k = 1, 2, 3) since in equilibrium and in the comoving frame Mq = 0. Here 
the naught subscript denotes equilibrium values and 5 prefix the fluctuations around 
them. The linearized set of equations, which the fluctuations should follow according 
to the linear regression assumption, is given by 

5h + n 59 = (21) 

4 K^o + Po) 8u v + — 5n v + -J-ST V - 2 V (5^) - (56 ,„ - \ (L TT 5T U + F v )=0 

(22) 
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nCJT + ^59- [h$ (LttT' u + F v )\ . k = (23) 

To arrive at this set of equations the local equilibrium assuption has been used to 
express e = e(n,T) and p = p(n,T). C n is the specific heat at constant particle 
number density. We conveniently introduced a vector quantity F u as a wildcard 
which is associated with the hydrodynamic acceleration in Eckart's fomalism, with 
a gradient of number density in the one proposed in Ref. |8|, or some other gradient 
following Israel's early worl 
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c 23| or the calculations shown in relativistic kinetic 

22j. As was emphatically pointed out in previous work 



theory standard books 
i 

6j, the exponential growth of perturbations can be picked up simply by performing 
a Fourier-Laplace transform followed by the calculation of the curl of Eq. ( 1221) . This 
isolates the transverse mode and allows for the assesement of its behavior in time. 
If F v is associated with a time derivative in the velocity, the procedure yields a 
cuadratic equation for the Laplace variable with one postive root which in turn leads 
to the unphysical growth of fluctuations. However, if F u is indeed a gradient of 
a thermodynamic quantity, as predicted by kinetic theory, its curl vanishes and a 
first order equation is obtanied. This yields an exponential decay for the transverse 
fluctuations with which one concludes that the instability of the equilibrium state is 
not longer present. For the longitudinal mode one obtains 



p Se + — V 2 5n + ^-V 2 5T 

—AV 2 59 — — ^V 2 <5T — \f v v =0 (24) 

Equations (|2ip . (123]) and f l24"|) constitute a set of three coupled equations for 5T, 
5n and 59. Taking their Fourier-Laplace transform a set of algebraic equations is 
obtained which may be solved to yield a dispersion relation from which its roots may 
be obtained. After a rather long calculation by inverse transformation one obtains 
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explicit solutions for 5n(f,t) and ST(r,t). 

As it is well known the normalized autocorrelation function of the density fluctu- 
ations is proportional to the dynamic structure factor 
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1* 



This calculation, 



in the non-relativistic case, yields the well known Rayleigh-Brillouin spectrum pre- 
dicted by Landau and Placzek in 1934 1 301 and measured experimentally for Argon 

I I 

in 1966 by Boon et al. [17| . This is one of the many experimental confirmations of 
the linear regression assumption. 

What we expect to find in the relativistic case is precisely a relativistic correction 



28], the 



to this spectrum. This does not occur when if oc w M . As shown in Ref. 
spectrum simply does nor exist. On the other hand, if F 11 cx Wn )V we do indeed 
recover the spectrum with relativistic corrections to Rayleigh's peak, which in the 
non-relativistic limit reduces to the classical form. Even if the experiment is not 
or cannot be performed for technological reasons, the result speaks in favor of a 
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23] and not Eckart's one. 



constitutive equation in the form given in Refs. js), Q, 
It is precisely due to the structure of Eq. (1T51) that one can show that the resulting 
linearized relativistic hydrodynamic equations are in full agreement with Onsager's 
linear regression assumption . This in turn implies that the equilibrium state 

is thermodynamically stable. 

V. DISCUSSION 

In the previous sections, we have reviewed recent results in relativistic LIT and 
kinetic theory and pointed out some key theoretical arguments related to them. In 
Sect. II, we emphasized the meaning of the generic instabilities found in Ref. Q 
as violations of Onsager's linear regression assumption. This hypothesis, as dicussed 
above, lies deep inside the theoretical setup of LIT and its violation is reason enough 
to discard a theory of irreversible processes. Based on this violation by Eckart's 
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linear theory, higher order formalisms were 



developed and now generally favored 



31 



However, based in our recent 



even though they also have been objected 
results and as outlined here in Sects. Ill and IV, the source of this unphysical 
behavior is the coupling of heat with acceleration assumed by Eckart's theory which 
is in disagreement with standard kinetic theory. These facts do not rule out the 
possibility of considering higher order, or extended, theories 33M35] but questions 



the need for them. That is, the unphysical behavior found in Eckart's formalism is 
corrected once kinetic theory is used in order to construct a heat flux tensor. These 
findings reopen the possibility of modeling relativistic fluids with first order theories. 

In Sect. II we outlined the calculation of Ref. js| where we obtained a constitutive 
equation that couples the heat flux with the thermodynamic forces VT and Vn. This 
is clearly consistent with the generalization of Onsager's theory since it couples fluxes 
with thermodynamic forces of the same tensorial rank, the former in this case being 
the heat flux and the latter WT^ and Wn^ in the (n, T, u^) representation. This 
constitutive equation is then introduced in Sect. IV to argue that with it no stability 
nor causality issues arise in the relativistic transport equations. The details of these 



28 



We also comment on the causality 



36|. 



calculations can be found in Refs. [6] and 
and refer the reader to a recent publication in this subject 

This work also serves to update a discussion that has been going on for some time. 
The first issue concerns the nature of the time component of heat flux four-vector. 
It is important to notice that we are using a projector in the constitutive equation 
(I15jl . This was questioned in previous work {25 1 but has since been analyzed. Indeed, 



in Ref. 



21] we show that a generalized relativistic Chapman-Enskog procedure for 



the relativistic Boltzmann equation is possible and consistent with irreversible ther- 
modynamics. We concluded that the fourth component of the heat flux as predicted 
by such a formal approach is still the internal energy itself. Thus, no "disspation in 
the time direction" is present in the theory and the projector should be included in 
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the constitutive equation as done by most authors. 

Secondly, the presence of heat in the stress-energy tensor must be addressed here 
following the discussion in Refs. J26 1 and J27 1. In Sect. IV, in particular in Eq. (1201) . 
we clearly state that here as in some of our recent work, the stress-energy tensor con- 
sidered is indeed the one proposed by Eckart strictly from a mathematical argument. 
However it still lacks a kinetic justification in the laboratory frame, eventhough some 



work along this line can be found in Ref. 



37fl. Nevertheless the relativistic heat terms 



can be readily shown to be present in the tensor in the comoving frame where most 
of our calculations are performed. Indeed we have lately adopted such a tensor as 
the fluid's stress energy tensor where, as mentioned earlier, the heat flux is included 
and the first term corresponds to an internal energy flux. Both energy and momen- 
tum balances are obtained from the general conservation law T.^ u = as pointed out 



above. On the other hand, the work in Ref. 



25] . which raised the discussion leading 



to Refs. [26] and 27], used a Meixner-like approach where the heat is not included 
in the stress energy tensor but is introduced jn the system via the construction of a 
total energy flux. The predictions found in 



28] finally lead to the criteria which, on 
experimental grounds, should prevail to decide which of the approaches is the correct 
one : Eckart's [5|, the one proposed by Sandoval- Villalbazo and Garcia-Colin in Ref. 



25] . or the one in Ref. |8j. This is still an open question. 



n 

It has to be cl arify ed here that in the present work as well as in Refs. [8j and 
Section IV of Ref. (28], the tensor given in Eq. ( 120]) is used as the fluid's stress energy 
tensor eventhough a kinetic proof of its structure, other than in the commoving frame, 
is still to be established. We consider that the question of whether this same tensor 
is the one to be included in Einstein's field equation is still open and should be 
addressed. The actual solution of such equation including heat sources should be 
critically analyzed in order to asses the possible effects of including such terms. This 
is clearly outside the scope of this review and will be addressed elsewhere. 
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